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Abstract 

We show that the potential wells V{x) = + ax^ + bx^ + cx"^ + dx"^ + f /x"^ with 
a central spike possess arbitrary finite multiplets of elementary exact bound states. 
Their strong asymptotic growth implies an ambiguity in their VT symmetrically 
generalized quantization (via complex boundary conditions) but the three eligible 
recipes coincide at our exceptional solutions. 
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1 Introduction 



Realistic calculations in quantum physics and field theory are often guided by a 
parallel study of a simplified quantum-mechanical model in one dimension. A typical 
example may be found, e.g., in ref. |1| where a self-interacting scalar field theory is 
being modified in non-perturbative manner. It is underlined there that one has to 
pay due attention not only to the interaction potential V{x) itself but also to the 
boundary conditions imposed upon solutions in infinity. 

In general the latter conditions are not unique. One of the most transparent 
explicit examples of their ambiguity has been offered by Bender and Turbiner p[. 
In essence, they considered a partially solvable sextic potential V{x) 
and its zero-energy ground-state wave function il){x) = exp(— a;^/4). In this model 
it is obvious that once you start working in the whole complex plane, x = ge'^ 
the asymptotic normalizability of the wave function keeps satisfied not only on the 
real line but, in general, within the four different asymptotic wedges defined by 
the elementary equation Re {—x^) < 0. This gives the admissible angles if G Sk, 
/c = 1, 2, 3, 4 lying within the four separate intervals. 



Schematically, the situation is depicted in Figure 1. One imagines that the eligible 
complex physical coordinates can be arbitrary curves with ends which do not enter 
the "forbidden" asymptotic domains. Thus, the current real line starts at ipi ^ S3 in 
the left infinity and ends at (fn G Si. Its slightly bizarre "Wick-rotated" alternative 
5*4 + 5*2 has also been discussed in the literature as a weird example of a system 
whose real spectrum is not bounded below The elementary zero-energy bound 
state itself is shared by both these non-equivalent spectra. 

From a retrograde point of view the choice of the sextic V{x) proved unfortunate 
since its menu (|ip is not yet rich enough. Further progress has only been achieved five 



Si = (-7r/8,7r/8), 
^3 = (77r/8,9 7r/8) 



^2 = (3 7r/8,5 7r/8), 
^4 = (ll7r/8,13 7r/8) . 



(1) 
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years later when Bender and Boettcher [^] came to the conclusion that a privileged 
role must be played by the pairs of sectors with a mirror left-right symmetry. This 
symmetry keeps the trace of its origin in field theory and is called VT symmetry. 
In the quantum mechanical context the Hamiltonians H have to commute with the 
product of V (parity) and T (complex conjugation or time reflection) |^. It is 
currently believed that this guarantees the reality of the spectrum for many non- 
Hermitian complex Hamiltonians [^]. 

In the very special class of these models V{x) = {ixY^ the complex plane is 
to be cut upwards in order to keep the picture unique. Bender and Boettcher |^ 
started from the smallest exponents 5 and picked up the lower sectors 

Sl = (-3 a - 7r/2, -A - 7r/2), 5^ = (A - 7r/2, 3 A - 7r/2) 

with the half- width A = 7r/(4 + 2 5). Of course, for any 5 > 1/2 there emerges an 
alternative pair of sectors 

Sl2 = (-5 a - 7r/2, -3 A - 7r/2), Sr2 = (3 A - 7r/2, 5 A - 7r/2) 

which has been used by Buslaev and Grecchi in their study of asymptotically quartic 
potentials [0]. In general the latter possibility remains compatible with the current 
real coordinates in the interval of 5 G (1,3). 

For all the asymptotically power-law models the ambiguity of quantization is an 
interesting phenomenon. Numerically, this has been documented by several studies 
which made use of the limiting transition 5 ^ oo or of the absence of the cut at 
(5 = 1,2,... 1^. An even simpler form of the 5 dependence takes place at the positive 
integers Z = 1 + 6/2 in the potentials V{x) = x'^^^^ + 0{x'^^~^) since, asymptotically, 
their wave functions il){x) ~ exp(— x^'^/2Z) are symmetric on the real line. 

There appear the two new VT symmetric pairs of sectors {Sl, Sr) at each odd 
value of the integer Z. Thus, the 5 = 4 and Z = 3 decadic oscillator is the first 
model with an ambiguity of this type. This is illustrated in Figure 2 where the single 
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left-right pair of the asymptotic boundary conditions of Figure 1 is replaced by the 
triple choice. All the three angles ipi G Sij and ipn G S^j with j = 1, 2 or 3 are 
equally compatible with the normalizability of the exponential ~ exp(— a;®/6). 
With this motivation we shall consider here all the decadic polynomial potentials 
complemented by a central spike, 

V{r) = + ar^ + br^ + cr^ + cir^ + //r^ . (2) 

These forces contain as many as five independent coupling constants and their real 



forms have been studied by several authors in the literature Here, we shall 

re-consider these "spiked decadic" interactions within the generalized quantum me- 
chanics of Bender et al [Q. It in effect weakens the current Hermiticity of the Hamil- 
tonian to its mere VT symmetry. For this reason one can construct more solutions 
in principle. We shall see below that such an expectation is well founded, indeed. 

2 Decadic oscillators 

The results of study of Hermitian Hamiltonians with interactions were summa- 
rized by Ushveridze in sec. 2.4 of his monograph [|ll| on the partially (so called 



quasi-exact ly) solvable models. This summary implies that the Hermitian decadic 
force lies somewhere in between the numerous purely numerical models (for which 



"it seems absolutely unrealistic to find an exact solution" [jTT|) and quasi-exactly 
solvable models in the narrower sense (where one requires the algebraic solvability 
for a multiplet of K states at any pre- determined integer K). 

The latter category (say, type-I QES) lies already quite close to the harmonic 
and other completely solvable models. Its properties are best exemplified by the 
standard Hermitian sextic model (cf. eq. (1.4.13) and related discussion in ref. [0). 
The former extreme without any solvability is usually illustrated by the Hermitian 
quartic oscillator (cf. section 1.3 in ref. [11]). In this comparison, the "intermediate" 
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decadic models (e.g., eq. (2.4.5) in the Ushveridze's book) admit merely a few (K) 
exact bound states with a strongly limited multiplicity K < 5 Such a property 
(let's call it QES of type II) is entirely trivial at i^' = 1 and still quite easily achieved 
at the first iew K > 1 |g . 



2.1 VT symmetric regularization 

The best illustration of influence of the replacement of Hermiticity by the mere VT 
symmetry of the Hamiltonian is provided by the popular quartic oscillators. In 



ref. [|T3| these potentials were shown to belong to the type I of the QES category. 
Obviously, VT symmetry plays a crucial role in such an improvement of solvability. 

The picture will be completed in what follows. We are going to demonstrate that 
a VT symmetrization of eq. (0) leads still to the maximal, type-I form of the QES 
property. We shall see that for the complexified decadic oscillators the integer K may 
be chosen arbitrarily large. It is worth noting that the same enhanced solvability 
admitting the arbitrarily large multiplets remains freely applicable in any spatial 
dimension D. 

In a preparatory step of the explicit constructions let us remind the reader that 
in any ( = i-th) partial-wave projection the dimensional differential Schrodinger 
equation with a virtually arbitrary regular central potential reads 

_ rf2 ^ L(L + 1) ^ ^1 ^ ^ ^ = £ + p _ 3)/2. (3) 

In a more general perspective, we can add a centrifugal-like spike to the regular force. 
This is known to preserve the exact solvability of the spiked harmonic oscillator . 
In a close parallel, we do not need to change the notation too much, re-defining only 
the angular momenta L = L{f) in such a way that 

i(i+l) = /+f, + ^U,+ ^V (4) 
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As a consequence, equation (|^) may be assigned a pair of independent solutions with 
the well known behaviour near the central singularity ~ 

In such a setting the "forgotten possibility" of a suitable ansatz lies in the simulta- 
neous use of both these solutions in the (complex) vicinity of the origin, 

^{r) ~ Cir-^ [1 + 0{r^)] + C2r^+^ [1 + 0{r^)] . (5) 

Similar idea is slightly counterintuitive but it has already been used in several papers 
on the quartic oscillators with 5 = 1. In 1993, Buslaev and Grecchi paved the way 
by the mathematically rigorous example of introduction of the complex coordinates 
for a problem of the present type. They achieved a regularization of the centrifugal 
spike by a constant imaginary shift of the real axis, 

r = r{x) = x — ie, a; G (—00,00), e > 0. (6) 

This makes the centrifugal spikes smooth and fully regular at all x, 

L{L + l) ^ L{L + l){x + ief ^ L{L + l) 
{x-ief (a;2 + £2)2 ^2 ^^K^)- 

A longer discussion of some consequences of this type of the VT symmetric reg- 
ularization has been provided by ref. [III]. There, both the even and odd wave 



functions of the exactly solvable harmonic oscillator and other models were assigned 
their separate analytic continuations to the complex x. In accord with expectations, 
the spectrum of the energies remained real. 

In the present paper we shall try to follow the same pattern and imagine that the 
two-term ansatz (^ may remain compatible with the single Taylor-series expansion 
of ip{r) (say, in the powers of r^) whenever the ratio r^+^ jr^^ is itself equal to an 
integer power of r^. In the other words, under a suitable convention L + 1 > — /, 
our key assumption will read L + 1/2 = M where M can only be a positive integer, 
M = l,2,.... 
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We shall see below that this type of constraint will lead to significant simplifica- 
tions of the solutions as well as to their easier interpretation. Indeed, the trivial choice 
of / = implies that M = i—l+D/2 = 1. This mimics the four- dimensional s— wave 
or two-dimensional p— wave situation since, in both cases, L = 1/2. Similarly we ar- 
rive at the alternative choice between the six-dimensional s— wave, four-dimensional 
p— wave or two-dimensional (i— wave Schrodinger equation at M = 2, etc. 



2.2 Correct asymptotics and recurrences 

In 1998, Bender and Boettcher [|l^ discovered the partial solvability of the VT 
symmetrized quartic polynomial oscillators. In the present language this means 
that they just employed the ansatz of the type (^) at the particular L = 0. This 
has been accompanied by the asymptotically bent choice of the complex integration 
path r —>■ r{x) G (p. In a way inspired by ref. |^ the latter construction has been 
extended to all L in our recent remark [ p!5|] . 

It is amusing to summarize that except the pioneering sextic study and its 



harmonic-oscillator simplification [HI WM, virtually all the available papers on the 



(partially) solvable VT symmetric polynomial potentials pay an exclusive attention 
to the asymmetric models V{r) ^ V{—r) of degree two [|16| and four [^, |13|, |15[. 
In this way all of them avoid the ambiguity problem but necessitate the complex 
couplings and acquire the counter- intuitive property Im V{t) ^ even on the real 
axis of coordinates t. 

The unpleasant asymmetries disappear within the decadic model @ which is 
such that V{r) = V{—r). Its asymptotic growth is also steeper than in the current 
solvable models [|l^. In what follows, we are going to show that this model represents 
in fact the "missing" last item in a list of all the quasi-exactly solvable polynomial 
models. The first step towards this not quite predictable result lies in the manifestly 



6 



normalizable ansatz 

/ ^6 4 2\ N-1 

Hr)=ew[-j-aj-0-j E /.^r^""^. (7) 

The insertion of this formula converts our differential Schrodinger equation (|^) into 
the finite set of recurrences 

Anhn+1 + Bnhn + Cnhn-1 + Dnhn-2 = 0, n = 0,l,...,N. (8) 

The use of the asymptotically optimal WKB-inspired parameters 

a = 2a, b = a'^ + 2p, c = c{N) = 2a(3 + 2M - AN - 2 

enables us to simplify the coefficients significantly, 

An = {2n + 2)(2n + 2~2M), Bn = E - f3 {An + 2 - 2M), 
Cn = - d- a{An-2M), D„ = 4 (A^ + 1 - n). 

This is the concise, linear algebraic formulation of our present problem. 

3 Terminating solutions 

Centrifugal spike in eq. binds the integers M to its strength / 7^ in a way 
prescribed by formula @. This means that certain non- vanishing spikes can emerge 
as the simple functions of the dimension and angular momentum, 

f = M^ - (£ - 1 + D/2)l 

In particular, the most elementary choice oi D = M = 1 and i = implies that we 
just have to fix / = 3/4. This value is, by the way, precisely equal to a boundary 
of a certain mathematical regularity domain (cf. our recent remark |]TB[ for more 
details). 
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3.1 Sturmian multiplets of couplings at M = 1 



We may notice that in accord with our definitions (^) two coefficients vanish com- 
pletely, A^i = and D^+i = 0. This is vital for the consistency of our ansatz (^. 
The third consequence of our assumptions reads Am-i = and has no obvious inter- 
pretation. Seemingly, it is redundant. Let us now pay more attention to its crucial 
and beneficial role. 

Starting from the first nontrivial choice of M = 1 we get simply Aq = 0. In 
such a case we can fix E = and discover that the whole over-determined set 
of our recurrences (H) degenerates to the mere square-matrix recipe. As long as 

= Cn{d) = (3"^ - a {An - 2M) - d we can write 



Ci(0) 5i 

D2 C2(0) 



Ai 



\ 



D 



N-l 



V 



ho 
hi 
h2 



\ 



h 



N-l 



( 



d- 



J 



ho 
hi 
h2 



\ 



h 



N-l 



(10) 



Bn-2 An-2 
Cn-i{^) Bn-1 
Dn Cjv(O) ) 

We see that a more or less routine diagonalization of a four-diagonal matrix deter- 
mines in principle the N different eigen-couplings d = dk with k = 1,2, . . . , N. These 
values may be found numerically at an arbitrary A^. In the other words, the whole 
problem becomes quasi-exactly solvable of type I. This is our first important result. 

It makes sense to introduce a shifted coupling F = d — (3'^ + 2N a. For the first few 
smallest dimensions > this simplifies the formulae and leads to the transparent 
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implicit definitions of tlie sliifted couplings F{d, N), 



F 







7V = 1, 



F2 + 16/3 -4^2 







= 2, 



-F^ + {16 a'^ - 64 /3)F + 256^0, N ^ 3, 
F^ + (160 /3 - 40 q;2) F^ - 1536 F + 144 - 1152 (3a'^ + 2304 /S^ = 0, AT = 4, 
-F^ + (80 q;2 - 320 /3) + 5376 + (8192 /3 - 1024 - 16384 P'^) F 
+196608 /3 - 49152 q;2 = 0, = 5, 

We may summarize that at M = 1, our solutions remain purely non-numerical up 
to the degree N — A. In a sufficiently broad part of the {a, (3)— plane (or, if you 
wish, of the (a, b)— plane of the octic and sextic coupling constants) the multiplets of 
exact and elementary bound states exist and are numbered by their (real) quadratic 
couplings dn at any value of the multiplicity N. 

3.2 Multiplets of energies at M = 2 

We have seen that our M — 1 multiplets have been formed by the so called Sturmian 
solutions. The role of their energy was marginal and fixed to the single value E — 0. 
Returning now back to our main story we have to move to the next integer M = 2. 
This replaces the above-mentioned choice of F = (i.e., in effect, of Bq = 0) by the 
virtually equally efficient simplification 



det 




which leads immediately to the compact formula 



d = d{E) 



F2 
T' 



M = 2. 
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The insertion of this energy- dependent harmonic strength d returns us back to the 
square-matrix secular equation ([To|) . With its three innovated diagonals 



An= {2n + 2){2n-2), Bn = E - (3 {An - 2), Cn = - E^A - a {An - 4) 

it defines the spectrum En in the very similar manner as above, i.e., as roots of a 
certain polynomial. We can display its first nontrivial sample, 

-E^ -2/3^4+ (-48a + 8/?2) + (-96/? a + 192 + 16/53)^2 
+ (-256 P - 512 + 192 p^a - 16 p^) E 
-1024 pa^- 1280 P^ + 4096 a - 32 p^ + 384 p^a = 0, 
M = 2, N = 3. 

Its roots are not non-numerical anymore but three of them remain manifestly real 



a.t a = P = where their form remains closed, (i?i, i?2, -E3) = (0,0, vT92). The 
consequent and precise specification of the whole domain of the reality of these roots 
remains as numerical a task as, say, its parallel studied in ref. [|T3| . 



3.3 More complicated multiplets at M = 3 etc. 

Explicit formulae with the integers M > 3 become appreciably more complicated. 
For all the really large truncations N ^ M they may remain useful. Their derivation 
would proceed along the same lines as before. One starts from the general pre- 
conditioning requirement 



det 



Bo Ao 
Ci B, 
D2 C2 



Ai 

Bm-2 Ai\i_2 

E>M-i Cm-1 Bm-i 



:iii 
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and its solutions have to be inserted in the, presumably, much larger main secular 
determinant 



det 



D2 C2 



V 



A, 








Bn-2 


An-2 


Dn-1 


Cn-1 


Bn-1 




Dn 


Cn 



. 



:i2) 



I 



The procedure only becomes inefficient at the larger integers M. In such a setting it 
would be more appropriate to treat both the "small" and "large" secular equations 
(pJ]) and (|T2|) on an equal footing, as a mutually coupled algebraic system. 

In the M ^ N setting one suddenly loses the main advantage of our present 
construction, viz., its reducibility to the single secular equation. In practice, one also 
has to replace the insertions oi E 01 d by the use of the so called Grobner bases. 
The algorithm determines both d and E at once and is routinely provided by the 



languages like MAPLE p 



Although the generalized M > 3 formulae and calculations need not necessarily 
become hopelessly complicated (and, in fact, give the nice results, e.g., in the M 00 
limit [^) their use would definitely require another motivation. Once we get that 
far, we would have no reason for maintaining our main assumption (1). At any real 
quantity M = 2L + 1 (reflecting a free variability of the coupling /) we would only 
have to replace the "small" dimension M in eq. ([Tl|) by the overall truncation N 
itself. 



4 Concluding remarks 

Results of our explicit constructions share their transparent algebraic character: 
Mathematically, one eliminates the auxiliary or "redundant" root of "the simple" 
condition ( pH] ) and ends up with the single "effective-Hamiltonian-like" square-matrix 
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eigenvalue problem (JT^)- This is the main merit (and general feature) of the QES 
systems of type I. Indeed, once we return to the explicit Hermitian decadic N = 2 
construction of sec. 2.4 in ref. we discover that, in the same language, the main 
shortcoming of the type-II QES constructions lies precisely in the necessity of solving 
several coupled eigenvalue problems at once. 

Our present sample non-numerical constructions available at the first few small- 
est integers M parallel in fact many other quasi-exactly solvable models. In par- 
ticular, the Sturmian constant-energy form of the elementary multiplets charac- 
terized already the very first quasi-exactly solvable (viz.. Coulomb plus harmonic) 
model as discovered by A. Hautot in 1972 |^. More complicated relation between 
the couplings and energies characterizes, e.g., the partially solvable anharmonicities 
~ {l + gr"^)"^ revealed by G. Flessas in 1981 p2| and explained algebraically in 1982 

In our present paper the most elementary multiplets using auxiliary integer M = 1 
proved purely Sturmianic. Their spike-shaped short-range part of the interaction ~ / 
is not arbitrary and can only vanish in even dimensions. Vice versa, the emergence of 
a spike for our s— wave multiplets in three dimensions is reminiscent of the so called 



conditionally solvable models where the choice of / = 3/4 is obligatory |T8 



At any auxiliary M , in comparison with all the Hermitian QES-I systems related 
to certain tridiagonal matrix representations of Lie algebras [^], all the decadic 
examples are distinguished by the four-diagonal and A^— dimensional secular deter- 
minants ([T2|). In the Hermitian setting this "four-diagonality" was in fact the main 
reason of the restricted type-II solvability. We have shown that only the appropriate 
complexification can move the decadic systems to a higher, type-I QES group. At 
present, this group which exhibits the appealing VT symmetry already encompasses 
the asymptotically decreasing quartic forces of Bender et al JT3| |TB[ and the Coulomb 
-|- harmonic model |jT6[ . 
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Figure captions 

Figure 1. Permitted domain for sextic oscillators 
Figure 2. Permitted domain for decadic oscillators 
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Figure 1 . Permitted domain for sextic oscillators 
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Figure 2 . Permitted domain for decadic oscillators 




